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ABSTRACT 

u. T ! ,e u X ^ ra r y/radi ° transient superluminal jet source GRO J1655-40 was recently suggested to contain a 
black hole from optical observations. Because it is a relatively close-by system (d ~ 3.2 kpc), it can likely 
provide us with rich information about the physics operating in both Galactic and extragalactic jet 
sources. We present the first simultaneous broadband high-energy observations of GRO J 1655— 40 
during the 1995 July-August outburst by three instruments: ASCA, WATCH /Granat, and BATSE/ 

r^r?\ .® nergy band fr ° m 1 keV to 2 MeV ° ur observations strengthen the interpretation that 

GKO J1655— 40 contains a black hole. We detected a two-component energy spectrum, commonly seen 
from other Galactic black hole binaries, but never detected from a neutron star system. Combining our 
results with the mass limits derived from optical radial velocity and orbital period measurements, we 
further constrain the mass of the central object to be between 3.3 and 5.8 Af e , above the well-established 
mass upper limit of 3.2 M Q for a neutron star (the optical mass function for GRO J1655— 40 is 
3.16 + 0.2 M q ). This system is therefore the first Galactic superluminal jet source for which there is 
strong evidence that the system contains a stellar mass black hole. The inclination angle of the binary 
system is constrained to be between 76° and 87°, consistent with estimates obtained from optical light 
curves and radio jet kinematics. 

Subject headings, binaries: close black hole physics — gamma rays: observations — 
stars: individual (GRO J1655-40) — X-rays: stars 


1. INTRODUCTION 

Radio observations are now beginning to show that rela- 
tivistic jets may be a more common feature of X-ray binary 
systems than previously thought. Some, such as SS 433 
(Margon 1988), Cygnus X-3 (Strom et al. 1989), and the 
more recently discovered superluminal transients GRS 
1915+105 (Mirabel & Rodriguez 1994) and GRO 
J1655-40 (Tingay et al. 1995; Hjellming & Rupen 1995) 
show proper motions implying velocities V ~ 0.2 0.9c. Of 
those, only GRO J1655— 40 has a compact object whose 
dynamically estimated mass is greater than 3 M 0 (Bailyn et 
al. 1995b). It is thus generally believed to be a black hole 
binary (BHB). Because it is a BHB with relativistic jet ejec- 
tion episodes, GRO J1655-40 is a good analog to the 
active galactic nuclei (AGNs), which are thought to contain 
supermassive black holes and frequently are observed to 
have relativistic jets associated with the central engine. 

The Galactic jet sources are of great interest because the 
dynamic timescales of BHBs are a factor of ~10 6 10 9 
smaller than those of typical AGNs (Shakura & Sunyaev 
1976). Therefore, dynamical variations of AGNs on time- 
scales from many years up to the Hubble time correspond 
in BHBs to intervals of seconds to years. Observations of jet 
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formations and decay in coincidence with studies of overall 
flux variations and spectral evolution allow the possibility 
of identifying the emission regions at various wavelengths. 

Studying the inner part of the accretion disk very close to 
the central black hole is critical for understanding the 
central engine and the jet ejection. This region can be 
observed at X-ray and hard X-ray wavelengths for BHBs. 
For AGNs, however, this region is much more difficult to 
observe. This is because the typical photon energy radiated 
from the inner accretion disk region of an AGN is much 
lower (~eV). Absorption of such low-energy photons in the 
complex surrounding environment and the interstellar 
medium is severe. Thus GRO J1655-40 is a unique system 
for understanding the physics involved in accreting black 
holes and relativistic jets in our own Galaxy and AGNs. 

Since the discovery of GRO J 1655 -40 (Zhang et al. 
1994a) by the Burst and Transient Source Experiment 
(BATSE) aboard the NASA Compton Gamma-Ray Obser- 
vatory ( CGRO ) on 1994 July 27, it has been extensively 
observed in radio, optical. X-ray, and gamma-ray bands. 
Radio observations have revealed several relativistic ejec- 
tion episodes (Tingay et al. 1995; Hjellming & Rupen 1995) 
from this system, following hard X-ray (above 20 keV) 
activity (Harmon et al. 1995). Optical observations (Bailyn 
et al. 1995a, 1995b) have determined the orbital period (2.62 
days), the nature of the companion star (F-type), and the 
mass limit of the central compact object ( > 3 M 0 ). 

The BATSE experiment detected four major outbursts 
from this system separated by about 120 days (Zhang et al. 
1995a). Two of them were also observed by the WATCH 
instrument on the Granat satellite (Sazonov, Sunyaev, & 
Lund 1996). The energy spectrum above 20 keV when the 
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source is in outburst is well described by a power law 
(photon index 1.5 3.1) up to at least 200 keV. OSSE /CGRO 
detected emission up to 600 keV without any detectable 
deviation from the power law (Kroeger et al. 1996). Searches 
in both BATSE and OSSE data for rapid variability 
resulted in upper limits of about 5% of integrated rms noise 
from 0.01 to 1 Hz (Kroeger et al. 1996;Crary et al. 1996). In 
light of the possible optical eclipses (Bailyn et al. 1995a; but 
cf. van der Hooft et al. 1996), both BATSE and OSSE data 
were searched for orbital eclipses and/or modulation. No 
indications of periodic behavior were found (Zhang et al. 
1997; Kroeger et al. 1996). 

The TTM instrument (2-27 keV) aboard the Russian Mir 
space station also observed this source on several occasions 
(1994 September and October and 1995 February) between 
the main hard X-ray outbursts detected by BATSE. An 
ultrasoft spectrum, compatible with an accretion disk 
origin, was observed during most of these observations 
(Alexandrovich et al. 1994; Alexandrovich, Borozdin, & 
Sunyaev 1995). Sometimes a hard X-ray tail was seen by 
TTM and HEXE (15-200 keV) also aboard the Mir station 
(Sunyaev et al. 1997). X-ray observations of GRO 
J 1655 — 40 in the 1-10 keV range with the Japanese ASCA 
satellite have revealed very unusual spectral characteristics. 
On two occasions, 1994 August 23 and September 27, when 
the source flux in the BATSE energy band (20-100 keV) was 
very low ( < 30-50 mCrab), ASCA detected an energy spec- 
trum with a significant high-energy cutoff above 4 5 keV. 
Absorption-line features were seen in the region of the Fe 
K-edge (Ueda et al. 1997). 

2. MULTI-INSTRUMENT OBSERVATIONS 

The outburst in 1995 July-August was simultaneously 
observed from 1 keV to 2 MeV by ASCA (1-10 keV), 
WATCH /Granat (8-20 keV), and BATSE/CGRO (20-2000 
keV) instruments. The X-ray telescope (Tanaka, Inoue, & 
Holt 1994) aboard the Japanese satellite ASCA uses focus- 
ing X-ray optics to concentrate X-ray photons in the 1-10 
keV range onto its detector planes, the solid state spectro- 
meter (SIS), and the gas imaging spectrometer (GIS). The 
SIS was not usable for this observation owing to the bright- 
ness of the source, which caused telemetry saturation and 
pulse pile-up. The GIS was usable for energy spectral mea- 
surements, though significant dead time was present. Dead 
time was corrected for in estimating the incident source flux. 
The ASCA observation occurred between 1995 August 
15.45 and 16.18 (UT), during the rise to the second peak of 
this outburst. Overall flux variations during the nearly 1 
day observation were less than 10%. 

The WATCH instrument (Lund 1986) aboard the 
Russian Granat spacecraft detected this outburst while the 
Granat Observatory was scanning the sky. In the scanning 
mode, the motors of the WATCH collimators do not rotate, 
and the spinning of the satellite is used to produce modula- 
tion patterns. A detailed analysis of the data was done to 
separate the contributions of GRO J 1655— 40 and the 
nearby source 4U 1700 — 37 to the modulation curve. 
Observations during the X-ray eclipses of 4U 1700-37 and 
simulations of the instruments’s imaging capability gave 
proof that the obtained light curve of GRO J 1655 -40 is 
correct. 

GRO J1655-40 has been continuously monitored by 
BATSE’s Large Area Detectors (LADs) (Fishman et al. 
1989) using Earth occultation analysis (Harmon et al. 1992) 


and Earth occultation transform imaging (Zhang et al. 
1993, 1994b, 1995b) techniques. The uncollimated nature of 
the LADs makes other active sources also detectable during 
these outbursts. The data are thus simultaneously fitted 
with these sources included to determine the net contribu- 
tions from GRO J 1 655 — 40. 

3. BROAD ENERGY BAND SPECTRUM 

The ASCA-BATSE measured spectrum is shown in 
Figure 1. The ASCA data cannot be fitted by standard 
single-component models. However, they are consistent 
with a two-component spectral model, consisting of a 
power-law component and a multicolor blackbody disk 
(MCBD) component (Shakura & Sunyaev 1973; see 
Mitsuda et al. 1984 and Makishima et al. 1986 for the 
MCBD model). The flux detected by WATCH in the 8 20 
keV range during this time was 620 + 50 mCrab (here the 
Crab flux was determined contemporaneously with the 
GRO J1655-40 observation). The WATCH detection is 
also plotted in the same figure. Extrapolation of the BATSE 
power-law spectrum into the WATCH energy range gives 
~240 mCrab. The WATCH flux thus independently sug- 
gests the presence of some additional soft component, 
which should contribute ~380 mCrab in the 8 20 keV flux. 
The ultrasoft component observed with ASCA fitted by the 
MCBD model would produce ~400 mCrab in the 
WATCH energy band; therefore, the value of the flux mea- 
sured with WATCH is in a good agreement with a model 
consisting of a power-law component and a MCBD com- 
ponent. No other spectral component is needed to account 
for the 1-200 keV overall spectrum, except for a weak 
absorption feature in the Fe K-fluorescence energy band of 
the ASCA data (Ueda et al. 1997). The ASCA continuum 
model parameters for this two-component spectrum are as 
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Fig. 1.— Joint ASC A/BATSE spectrum with WATCH integrated flux 
for 1995 August 15-16. Spectral models and data points are shown ( ASCA 
error bars are extremely small). The total ASCA energy spectrum is com- 
posed of two components, i.e., an ultrasoft component described by a 
multicolor blackbody disk (MCBD) model and a power-law component 
that is consistent, when extrapolated into the BATSE energy range, in both 
photon index and flux with the independently determined BATSE power 
law. The WATCH data point is consistent with the sum of these two 
components. No third spectral component is required for the entire 1-100 
keV spectrum, except for a weak absorption-line feature in the iron K- 
fluorescence energy band of the ASCA data. 
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companion mass of 0.1 M 0 (below which such a binary 
system would not be formed), a sufficiently high inclination 
angle would result in the center of the accretion disk to be 
eclipsed by the companion. This “eclipse” line is depicted 
by the nearly vertical thin solid line. The parameter space on 
the right side of this line is a part of the eclipsing zone . For a 
lower inclination angle, a more massive companion will 
also result in such eclipse. The “eclipse” line corresponding 
to 0.32 A/ 0 is shown as the vertical thin dotted line. The 0.32 
M 0 is chosen because this line intersects with the optical 
M-i relationship for the same companion mass on the upper 
bound of the X-ray M-i. Therefore, if the system parameters 
fall in the area B, the accretion disk center will be eclipsed 
by the companion. In area A, no such eclipse will happen. 
The mass of the companion cannot be higher than 0.32 M 0 
for a noneclipsing system. 

A reasonable mass of the companion is between 0.5 and 
1.5 M 0 . In this case, this system will be the first detected 
eclipsing BHB in our Galaxy. Eclipse mapping may be 
applied to determine the details of the inner disk region for 
the first time. The BH mass lies between 4.1 and 5.3 Af G , 
and the inclination angle is between 80° and 86?5. The inner 
disk radius can be estimated to be between 25 and 36 km, 
considerably larger than that of neutron star binaries ( ~ 10 
km). On the other hand, this system has an unusually large 
space velocity, which indicates that it probably has experi- 
enced an unusual evolution (Brandt, Podsiadlowski, & 
Sigurdsson 1995). Therefore, the companion could possibly 
be severely undermassive. Constraints from evolutionary 
theories suggest that the companion cannot be less massive 
than 0.23 M G (Brandt et al. 1995). From Figure 4, it is clear 
that the allowed parameter space between 0.23 and 0.5 M Q 
companion is very small. This system may also be an impor- 
tant object for our understanding of the formation and 
evolution of BHBs. 

Further support that the system contains a black hole is 
the lack of an additional thermal component commonly 
seen from a neutron star surface with a higher temperature 
(typically a ~ 2 keV blackbody spectrum, which is consider- 
ably harder than the MCBD spectrum, i.e., the ultrasoft 
component). The ultrasoft component by itself, however, 
cannot be taken as a signature of a black hole, since 4U 
0142 + 614, a neutron star pulsar system, has displayed such 
an ultrasoft spectrum (Mereghetti, Stella, & De Nile 1993; 
Israel, Mereghetti, & Stella 1993). Because of the mass acc- 
retion rate dependency of the disk temperature, a neutron 
star system may be observed with a very low disk tem- 
perature when its mass accretion rate is low and therefore 
the X-ray luminosity is low. We therefore argue that a lumi- 
nous (~ 1 L Eddington /or 1 M 0 ) ultrasoft X-ray spectrum is a 
very strong indication of a black hole system and should be 
used for selecting BHCs. A prominent power-law com- 
ponent associated with such a bright ultrasoft component, 
however, has never been observed from a neutron star 
system but is now seen from persistent, transient, and jet- 
type BHBs. Such a two-component spectrum may , therefore , 
be a firm signature from a black hole system (Sunyaev et al. 
1988, 1994; van der Klis 1994; van der Klis & van Paradijs 
1994; White 1993; Tanaka & Lewin 1995; Zhang et al. 
1996). 

6.2. Origin of the Power-Law Hard X-Ray Tail 

We are reasonably confident that the ultrasoft com- 
ponent originates from the accretion disk owing to the 


gravitational energy release of the accreted material in the 
X-ray energy band. The origin of the hard X-ray power-law 
component is, however, more controversial, although 
inverse-Compton upscattering of low-energy photons by 
fast-moving particles (electrons or protons) is commonly 
accepted as the radiation mechanism. The nondetection of 
short-term variations from the hard X-ray power-law com- 
ponent suggests a different origin of this power law than for 
the low-state hard X-ray power laws (with cutoff) observed 
from both black hole and neutron star binaries. Thermal 
Comptonization models normally used to interpret the low- 
state hard X-ray power laws cannot be applied here, regard- 
less of the nature and geometries of the hot Compton cloud, 
because the relatively steep power law extended to a very 
high energy implies a very high electron temperature (50- 
200 keV) and a very small optical depth (0.01 0.1). It is very 
hard to keep this optically thin cloud sufficiently stable to 
account for the low (or no) short-term variability in the 
hard X-ray power-law component. Cooling of the high- 
energy electrons by the copious soft X-ray photons in the 
high state seems to be an additional complication, unless 
there is a powerful underlying heating mechanism. 

Based on the detected relativistic jets from this source 
and the correlation of radio flares with hard X-ray activity, 
it has been suggested that the fast-moving particles 
responsible for the radio emission from the jets (Levinson & 
Blandford 1996) also scatter low-energy photons via the 
inverse-Compton process to produce the observed hard 
X-rays. We believe this is unlikely since this and the pre- 
vious outburst were not detected to have any radio emis- 
sion, unlike the first two outbursts with strong radio 
emission. The overall radio emission decayed exponentially 
since the initial hard X-ray outburst, while the hard X-ray 
outbursts with similar peak luminosities occurred 4 times 
since the initial outburst with a separation of about 120 
days. Absence of the jet ejection and any radio emission at 
all during the last two hard X-ray outbursts argues strongly 
against the scenario in which hard X-rays are produced 
directly or indirectly from the jets. 

The upper limit of about 0.5 mJy at 3.6 cm (Hjellming 
1996) allows us to set a limit of the size of the possible radio 
emission region. The surface brightness (due to incoherent 
synchrotron radiation) temperature upper limit of 10 12 K 
(Shu 1991) indicates that the dimension of the possible radio 
emission region is less than 1.7 x 10 1 1 cm (R. M. Hjellming, 
private communication). This is smaller than the size of the 
binary system (~ 10 12 cm) but can still be several orders of 
magnitude larger than the soft X-ray (blackbody) emission 
region (~ 10 7 cm). Since this region is so close to and very 
likely covers completely the soft photon emission region, 
the high-energy electrons required for producing synchro- 
tron radiations will be cooled down effectively by the 
inverse-Compton scattering process, thus preventing signifi- 
cant synchrotron radiations. Therefore, the power-law hard 
X-ray photons cannot be produced via the synchrotron 
radiation mechanism. 

Another model involves the converging flow of fast- 
moving material near the black hole horizon (Chakrabarti 
& Titarchuk 1995). This model does not have the difficulties 
discussed above for other models. In this model, the power- 
law component is produced via the inverse-Compton scat- 
tering of low-energy photons of the ultrasoft component by 
the bulk motion of the converging flow. This model predicts 
a spectral break energy of ^511 keV divided by the total 
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mass accretion rate in terms of the Eddington rate corre- 
sponding to the central compact object mass. The detected 
luminosity during the ASCA observation is ~ 10 38 ergs s“ 1 
and could be a factor of 2-5 higher at other times (scaled up 
from the detected WATCH and BATSE fluxes). Therefore, 
the Eddington accretion rate is <1 for a 3-5 M 0 black 
hole. The unbroken power law up to 200-600 keV is cer- 
tainly consistent with this prediction. It is also interesting to 
note that a relativistic converging flow cannot take place 
near a neutron star surface when the mass accretion is high 
owing to the radiation pressure from the neutron star 
surface. This provides a possible explanation why such a two- 
component spectrum has never been observed from a neutron 
star binary . We consider, however, that it is still difficult 
within this model to explain the range of spectral variations 
during this outburst, especially the very hard power law 
with a huge thermal excess at low energies near the end of 
this outburst. If GRO J 1655 -40 were indeed an eclipsing 
system, this model would be excluded, since it requires the 
converging flow and thus the hard X-ray production region 
to be very close to the BH. However, no hard X-ray eclipse 
or any significant hard X-ray orbital modulation is seen in 
the BATSE or OSSE data. 

6.3. Origin of the Short-Term Intensity Variations 

Short-term variability was detected for the first time from 
a Galactic superluminal jet source in a high state. In view of 
the low (or no) variability in the 20-100 keV band, there are 
at least two possibilities to account for the detected short- 
term variability in the 1-10 keV band. One is that the ultra- 
soft component is variable. We consider this to be unlikely 
since no similar variability is seen from the ultrasoft com- 
ponent of other BHBs, unless the nature of this ultrasoft 
component is different owing to the jet ejection nature of 
this system. (However, the energy spectrum does not look 
different from other BHBs.) Two other BHCs, namely GX 
339-4 (Miyamoto et al. 1991) and GS 1 124 — 68 (Miyamoto 
et al. 1991, 1994; Ebisawa et al. 1994), in their high state 
have been detected with a combination of strong hard 
X-ray power-law tails and moderate 5%) rms variations. 
These variations were found to be related to the power-law 
tails, and such a power-law like PDS from GRO J 1655 — 40 
was not seen from them. 

The other possibility is that the power-law component is 
variable but that the variability has an energy dependency 
so that the 1 — 10 keV power law is significantly variable but 
the 20-100 keV power law is not. To account for the 5.5% 
variability with only about 25% of the 1-10 keV flux in the 
power-law component, the intrinsic variability in the 1 10 
keV power-law component has to be as high as 20%. 
Energy-dependent variations have been observed pre- 
viously from other BHBs (Ebisawa et al. 1994; Miyamoto et 
al. 1991). The combined 1-10 keV variability and 20 100 
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keV low (or no) variability is, however, unusual among 
BHBs. 

Between the hard X-ray outbursts, almost identical short- 
term variability has also been observed by ASCA from 
GRO J 1655 -40 and GRS 1915+105 (Ebisawa et al. 1997). 
GRS 1915+105 is the only other detected Galactic super- 
luminal jet source. Its energy spectrum (power-law with 
cutoff above 5 keV, plus Fe-K absorption feature) is signifi- 
cantly different from the two-component spectrum we pre- 
sented here. This strongly suggests that this kind of 
short-term variability is related not to a specific energy spec- 
tral component but to the property that these systems can 
produce superluminal jets . This variability is, however, not 
directly related to the presence of jets because of the lack of 
significant radio emission during these observations. 

7. SUMMARY 

The broadband high-energy observations presented here 
have for the first time detected a two-component energy 
spectrum and short-term variability in the 110 keV band 
from a superluminal jet source. The former is shared by all 
types of identified BHBs but not by any known neutron star 
system, and it is therefore a firm signature from BHBs. The 
latter is unusual among BHBs and is probably somehow 
related to the uniqueness of this source as a superluminal jet 
BHB. We have constrained the mass of the central compact 
object and the binary system inclination angle using the 
X-ray energy spectrum combined with optical radial veloc- 
ity measurements. Our results further support the system 
containing a black hole with the highest inclination angle of 
all known BHBs in our Galaxy. The origin of the hard 
X-ray power-law component is still uncertain. Thermal 
Comptonization is unlikely to be the hard X-ray pro- 
duction mechanism, in view of the rather steep and 
unbroken power law up to very high energies. Jet origin 
models also seem difficult owing to the apparent lack of 
significant radio emission during this and the previous hard 
X-ray outburst. The converging flow model seems consis- 
tent with the spectral data, but it requires the hard X-ray 
production region to be compact. This can be tested if this 
system is an X-ray eclipsing BHB. The origin of the short- 
term variability in the 1-10 keV band is also puzzling. This 
is perhaps related to the unique nature of this source as a 
superluminal jet BHB. 
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Note added in proof . — After the acceptance of this article for publication, we became aware of new results obtained by C. D. 
Orosz and J. Bailyn (ApJ, in press [1997]) from optical light curve modeling of GRO J 1655 — 40. Their mass value 
(7.02 ± 0.22 M 0 ) and the binary inclination angle (69?50 + 0?08) are significantly different from what we obtained by 
combining the X-ray energy spectrum and the optical mass function. The apparent discrepancy might be caused by any of 
several factors. For example, if the black hole is rotating in the same direction as the accretion disk, the curves in our Figure 4 
would all move upward, which implies a higher black hole mass for a given inclination angle. Applying the inferred black hole 
mass value and the inclination angle from Orosz and Bailyn to our X-ray M-f relationship and assuming the black hole has 
nonzero angular momentum (Kerr black hole), we obtain an inner disk radius of ~2GM/c 2 , to be compared with 6GMjc 2 for 
a Schwarszchild ( nonrotating ) black hole. This inner disk radius would correspond to the last stable orbit of a ~95% 
maximally rotating black hole. Also, the X-ray spectral model we used may not be sufficiently accurate when the inclination 
angle is very high and the inner disk boundary is very close to the horizon of the black hole. 




